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ABSTRACT 

As revealed previously, under the assumption that some pulses of gamma- 
ray bursts are produced by shocks in spherical fireballs or uniform jets of large 
opening angles, there exists a standard decay form of the profile of pulses 
arising from very narrow or suddenly dimming local (or intrinsic) pulses due 
to the relativistic curvature effect (the Doppler effect over the spherical shell 
surface). Profiles of pulses arising from other local pulses were previously 
found to possess a reverse S-feature deviation from the standard decay form. 
We show in this paper that, in addition to the standard decay form shown 
in Qin et al. (2004), there exists a marginal decay curve associated with a 
local 8 function pulse with a mono-color radiation. We employ the sample of 
Kocevski et al. (2003) to check this prediction and find that the phenomenon 
of the reverse S-feature is common, when compared with both the standard 
decay form and the marginal decay curve. We accordingly propose to take the 
marginal decay curve (whose function is simple) as a criteria to check if an 
observed pulse could be taken as a candidate suffered from the curvature effect. 
We introduce two quantities A\ and A2 to describe the mentioned deviations 
within and beyond the FWHM position of the decay phase, respectively. 
The values of A\ and A2 of pulses of the sample are calculated, and the result 



suggests that for most of these pulses their corresponding local pulses might 
contain a long decay time relative to the time scale of the curvature effect. 
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1 INTRODUCTION 

In spite of the temporal structure of gamma-ray bursts (GRBs) varying drastically, it is 
generally believed that some well-separated pulses represent the fundamental constituent of 
GRB time profiles (light curves) and appear as asymmetric pulses with a fast rise and an 
exponential decay (FRED) (see, e.g., Fishman et al. 1994). 

Due to the large energies and the short timescales involved, the observed gamma-ray 
pulses are believed to be produced in a relativistically expanding and collimated fireball. The 
observed FRED structure was interpreted by the curvature effect as the observed plasma 
moves relativistically towards us and appears to be locally isotropic (e.g., Fenimore et al. 
1996, Ryde & Petrosian 2002; Kocevski et al. 2003, hereafter Paper I). Several investigations 
on modeling pulse profiles have previously been made (e.g., Norris et al. 1996; Lee et al. 
2000a, 2000b; Ryde et al. 2000, 2002). Several flexible functions describing the profiles of 
individual pulses based on empirical relations were presented. As derived in details in Ryde 
et al. (2002), a FRED pulse can be well described by equation (22) or (28) there. Using 
this equation, they could characterize individual pulse shapes created purely by relativistic 
curvature effects in the context of the fireball model. 

Qin (2002) derived in details the flux intensity based on the model of highly symmetric 
expanding fireballs, where the Doppler effect of the expanding fireball surface is the key 
factor to be concerned. The formula is applicable to cases of relativistic, sub-relativistic, 
and non-relativistic motions as no terms are omitted in the corresponding derivation. With 
this formula, Qin (2003) studied how emission and absorbtion lines are affect by the effect. 
Recently, Qin et al. (2004, hereafter Paper II) rewrote this formula in terms of the integral of 
the local emission time, which is in some extent similar to that presented in Ryde & Petrosian 
(2002), where relation between the observed light curve and the local emission intensity is 
clearly illustrated. Based on this model, many characteristics of profiles of observed gamma- 
ray burst pulses could be explained. Profiles of FRED pulse light curves are mainly caused 
by the fireball radiating surface, where emissions are affected by different Doppler factors 
and boostings due to different angles to the line of sight, and they depend also on the width 



and structure of local pulses as well as rest frame radiation mechanisms. This allows us to 
explore how other factors such as the width of local pulses affect the profile of the light curve 
observed. 

Revealed in Paper II, there exists only a slowly decaying phase in the light curve as- 
sociated with a local S function pulse, for which no rising phase can be seen. For a local 
pulse with a certain width, the light curve observed would contain both the rising and de- 
caying parts. It was revealed that light curves arising from very narrow local pulses and 
those arising from suddenly dimming local pulses share the same form of profiles in their 
decaying phase, which was called a standard decay form (see Paper II). For a common local 
pulse for which the decaying time is not short enough, the profile of the decay portion of 
the resulting light curve would significantly deviate from the standard form. It is interesting 
that, the deviation could be characterized by the feature of a reverse "S" (see Paper II Fig. 
5). We wonder if this indeed holds for FRED pulse GRBs. If it holds, how can we describe 
this deviation? Motivated by this, we explore quantitatively in this paper the deviation of 
light curves of gamma-ray burst pulses from the so-called standard decay form. A sample of 
FRED pulse sources will be studied. 

The paper is organized as follows. In section 2, we analyse the deviation of light curve 
pulses associated with gamma-ray burst spherical fireballs or uniform jets, from the standard 
form and a marginal curve. In section 3, we examine the deviation deduced from a sample 
detected by the BATSE instrument on board the Compton Gamma Ray Observatory. Dis- 
cussion and conclusions are presented in the last section. 



2 THEORETICAL ANALYSIS 

As derived in details in Qin (2002) and Paper II, the expected count rate of a fireball within 
frequency interval [vi, 1/2] can be calculated by 

2nR3 jre „ /(re)(l + _ t + re)dTe ^ go^e)^ 

fJ(q-\ ~fl,mm * /-i \ 

' ' ~ hcD^{\- 0)2(1 + _£_ r )2 '"" ~ ' ' 

which is just equation (21) in Paper II. This formula was derived under the assumption that 
the fireball expands isotropically with a constant Lorentz factor T > 1 and the radiation 
is independent of direction. Present in the formula, tq is a dimensionless relative local time 
defined by tq = c(tg — t c )/R c , where tg is the emission time (in the observer frame), called 
local time, of photons emitted from the concerned differential surface dse of the fireball (9 



is the angle to the line of sight), t c is a constant which could be assigned to any values of to, 
and R c is the radius of the fireball measured at t e = t c . In equation (1), I{tq) represents the 
development of the intensity magnitude in the observer frame, and go,u(^o,e) describes the 
rest frame radiation mechanism, with z/ o e being the rest frame emission frequency which is 
related to the observation frequency v by the Doppler effect. Variable r used in the formula is 
a dimensionless relative time defined by r = [c(t — t c )—D + R c ]/R c , where D is the distance 
of the fireball to the observer, and t is the observation time. As analyzed in Qin (2002) 
and Paper II, the relative time r is confined by 1 — cos# min + (1 — (3 cos 9 mill )r e min < r < 
1 - cos# max + (1 - (3 cos 6> max )T0 imax , and the integral limits r eMn and fe :max are determined 
by f e , min = max{T0 )min , (r - 1 + cos 6» max )/ (1-/3 cos 6» max )} and r ejmax = min{r ejmax , (r - 1 + 
cos6» min )/(l - /3cos6> mm )}, respectively, where we assign 6> min < 9 < 6> max and T e>min < r e < 

A local S function pulse, I{to) = cIo8{jg — Tg } o)/R c (where Iq and r^o are constants), 
when inserting it into (1), would produce an observed light curve of equation (35) in Paper 
II, which is 

CM = 2nR ' ! ° S i w ~ - Co(r), (2) 
with 

r , n _ (l+^ ) 2 (l + re,o-r) 

Co{ ) = r3(i- m i + ^ )2 • (3) 

Note that the range of r, within which the radiation of the local 5 function pulse over the 
concerned area is observable, is 1 — cos# min + (1 — f3 cos^mnjr^o < r < 1 — cos^max + (1 — 
/ 3cos^ max )r 6 i i o (see Paper II). As a product of the local 5 function pulse, light curves C(r) 
and C (t) reflect nothing but the pure curvature effect. For r>l, the term (3r/(l — (3) in 
equation (3) could be written as (3r/(l — j3) ~ (t — t )/(R c /2T 2 c), where to is a constant. We 
find that R c /2T 2 c is exactly the time scale of the curvature effect (see equation [5] in Paper 
I). (One can check that, in terms of local time, this curvature effect time scale becomes 

Rc/C.) 

In the case of the local 5 function pulse, as shown in Paper II equation (37), z/ ,e and v are 
related by z/ ,6» = [l+/5r/(l — — /3)ri//(l+/3T0 iO ), from which one gets dvQ y e/^o,e = dv/is. 
Thus, when taking go,v{vo,e) as a 5 function (i.e., when considering a mono-color radiation) 
and when interval [^1,^2] is large enough, C(r) would become C$(t) (differing only by a 
factor). Thus, Co(r) represents the light curve arising from a local 5 function pulse and 



associated with a mono-color radiation. As the radiation concerned (the GRB spectrum) is 
not a mono-color one and must last an interval of time, we call function C (t) as a marginal 
decay curve. 

As mentioned above, a sudden dimming or a very narrow local pulse could produce 
a standard decay form of light curves (see Paper II). According to Fig. 5 of Paper II, the 
standard form could be represented by equation (2). (Note that, a local 5 function pulse is an 
extremely narrow local pulse, and it belongs to the class of suddenly dimming local pulses.) 
Shown in Preece et al. (1998, 2000), a = —1 and f3 = —2.25 are typical values of the lower 
and higher indexes of the Band function spectrum, deduced from most GRBs observed. We 
thus define C{t) associated with the rest frame Band function spectrum of a$ = — 1 and 
/3 = —2.25 as a standard decay form which was mentioned in Paper II previously. 

With these two curves, we are able to explore how light curves associated with different 
local pulse forms deviate from standard ones, and with the deviation we might be able to 
estimate some parameters of local pulses. 

Taking T 9 fi=0, 9 min =0, and 6 max = n/2, we get < r < 1, and equation (3) becomes 

Co{r) = vHi-mil&T)*- (4) 

For the sake of comparison, we normalize light curves of (2) and (4) in intensity and re-scale 
its variable, r, by r = or' + b, so that the peak count rate is located at r' = and the 
FWHM position of the decay portion is located at r' = 0.2 (see Paper II). 

Formula (1) suggests that, except the state of the fireball ( i.e., T, R c and D), light curves 
of sources depend only on I{tq) and go, u {vofl)- We assume in this paper the common empirical 
radiation form of GRBs as the rest frame radiation form, the so-called Band function (Band 
et al. 1993) that could well describe spectra of most sources (see, e.g., Schaefer et al. 1994; 
Ford et al. 1995; Preece et al. 1998, 2000), and adopt different forms of local pulses, I(tq) 
in equation (1), which will produce different light curves, to study the deviation from the 
standard forms. 

Let us consider two local pulses. The first is a local pulse with an exponential rise 

I(t 9 ) = I CXp(— (T0 )min <T 6 < T0 )max ), (5) 

and the second is a local pulse with an exponential decay 

I {r e ) = I exp(-— — Z^EBE) ( r(9min < Te ). (6) 



Figure 1. -Plots of the normalized and re-scaled light curve C(r') arising from the local exponential rise pulse (5) (the left 
panel), the local exponential decay pulse (6) (the middle panel), and three other local pulses (7) (dash line), (8) (dot line), and 
(9) (dash dot line) (the right panel). The solid and short dash lines represent the standard decay form and the marginal decay 
curve in the left and middle panels, respectively. A Band function rest frame radiation form with ao = —1 and /3o = —2.25, 
and the frequency range of 100 < ^/hi,p< 300, are adopted, and we take T = 100, cr=0.02 (dash line), 0.2 (dot line), 2 (dash 
dot line), and 20 (dash dot dot line), respectively. 

Following Paper II, we take r 9tmax = 10a + Tg jmin (in this case the interval between TQ >max 
and r^min would be large enough to make the rising part of the local pulse close to that 
of the exponential pulse) and TQ >min = 0. Light curves arising from these local pulses are 
normalized and re-scaled in the way mentioned above, which are shown in Fig. 1. 

We find from Fig. 1 that there is no significant deviation of the light curves arising from 
local exponential rise pulses from the standard decay form in the decay portion of the light 
curve, just as what illustrated in Paper II. However, there are a significant "positive" devi- 
ation (over the standard form) of the light curves associated with local exponential decay 
pulses from the standard decay form within the range of r' E [0.0,0.2], and a significant 
"negative" deviation in the range of r' > 0.2, which we call a reverse "S" deviation char- 
acteristic. In addition, one finds from Fig. 1 that there exists a "negative" deviation of the 
standard decay form from the marginal decay curve within the range of r' > 0.2, but within 
the range of r' E [0.0, 0.2] there is no deviation between the two. 

We define the positive and negative deviation areas from the standard decay form (and/or 
the marginal decay curve) as A 1 and A 2 , respectively. Relation between Ai and A 2 and those 
between these areas and a are presented in Fig. 2 (here A 2 is calculated within the range 
of t' E [0.2,2.0], as light curves of local exponential decay pulses overlap each other in the 
range of r' > 2). As shown in Fig. 2, logA\ is linearly correlated with logA 2 . In addition, we 
find that logAi and logA 2 increase linearly with logo within the range of o < 0.05. A linear 
analysis produces logA x = —0.993 + 0.873loga and logA 2 = —0.278 + 0.952/ogcr. However, 



when a being large enough (say a > 2), A\ and A 2 would not change with a (in other wrods, 
they are saturated). 

To check if local pulses with different rising curves but sharing the same decaying curve 
would lead to much different profiles of light curves in the decay phase, we consider three 
other local pulses. The first consists a power law rise and an exponential decay: 



Tnmin ^ (r eMn <r e < r e ,o) 



Hje)=I Q { ^o-re,^' v - - — — (7) 

exp (_I^M) (j efi <Te < T(? max ) 

The second is an exponential rise and an exponential decay local pulse: 

?(„) = I { eXp( ^ ) (T » ■ £ T " £ r "> . (8) 

exp(-^M) („,„ < r„) 

The third is a Gaussian rise and an exponential decay local pulse: 

j, , j^M-C-^f] (r e , nm <r e <r efi ) 

exp(-^M) (r 9ja < r e ) 

We take a = 2.0 for the decaying part of these three local pulses. For the first one we take 
li — 2 and FWHM = 2.0, for the second and third ones we adopt o\ = 2.0. 

Light curves arising from these local pulses are presented in Fig. 1 as well. One finds 
that the three light curves possess the same decaying profile. This suggests that the profile 
of the decaying part of the light curve is determined only by the decaying curve of the 
corresponding local pulses. 

Besides the local pulses discussed above, we also study other forms of local pulses, such 
as a power law rise and power law decay pulse as well as a Gaussian local pulse. Relations 
between a or FWHM and A 1 and A 2 for the light curves arising from these local pulses are 
presented in Fig. 3. The same conclusions obtained above hold for these cases. 

According to the analysis above, one finds that: a) light curves associated with local 
pulses without a decaying portion would bear the standard decay form in their decaying 
phase, which was concluded previously in Paper II; b) there would be a reverse "S" feature 
deviation of the light curves arising from local pulses containing a decaying portion, from the 
standard decay form, which could also be concluded from Fig. 3 of Paper II; c) the deviation 
concerned could quantitatively described by areas A\ and A 2 defined above, and these two 
quantities are linearly correlated with the width of the decaying curve of the local pulse 
when the latter is small enough (say, when o < 0.05 in the case of an exponential decaying 
local pulse); d) all curves associated with a continuum spectrum (including the curve of the 
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Figure 2. -Plots of logAi and logA2 versus logo- (the left panels), and logA2 versus logAi (the right panels) for a typical 
local pulse. In the two left panels, the open square and open circle represent A\ and A2, respectively, associated with a local 
pulse with an exponential rise and an exponential decay. Quantities A\ and A2 presented in the two upper panels represent the 
deviation of pulses from the standard decay form, and those in the two lower panels describe the deviation from the marginal 
decay curve. Other parameters are the same as those adopted in Fig. 1. 
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Figure 3. —Plots of logA\ versus logFW H M or log a (the left panel), and logA2 versus logFW H M or logcr (the right panel) 
for some local pulses. Quantities A\ and A2 represent the deviation of the light curves associated with the local exponential 
decay pulse (the open square), the local power law decay pulse (with the power law index /i=2) (the cross), and the local 
Gaussian pulse (the open circle), from the standard decay form. Other parameters are the same as those adopted in Fig. 1. 



standard decay form) are well below the marginal decay curve beyond the FWHM position 
of the decay phase (say, within the range of t' > 0.2); e) within the FWHM position of the 
decay phase (say, in the range of r' G [0.0, 0.2]), the standard decay form and the marginal 
decay curve are not distinguishable. 



3 DECAYING FORM SEEN IN A FRED PULSE SAMPLE 



At least two questions urge us to employ a FRED pulse sample to make the following 
analysis. One is whether or not the reverse S-feature deviation characteristic indeed exists 
in the observed light curves, and the other is if the profiles of these light curves are indeed 
well below the marginal decay curve beyond the FWHM position of the decay phase. 

To study this issue, we utilize the light curves of the sample presented in Paper I, where 
the data are provided by the BATSE instrument on board the CGRO (Compton Gamma 
Ray Observatory) spacecraft. Of this sample, we find only the data of 75 individual pulses, 
which are employed in the following. It is well-known that pulses of a GRB show a tendency 
to self-similarity for each energy band (see, e.g. Norris et al. 1996). We thus consider in 
this paper only the light curve of channel 3 of BATSE, as signals in this channel are always 
significant. 

The background of light curves is fitted by a polynomial expression using 1.024s resolution 
data that are available from 10 minutes before the trigger to several minutes after the burst. 
The data along with the background fit coefficients can be obtained from the CGRO Science 
Support Center (CGROSSC) at NASA Goddard Space Flight Center through its public 
archives. All of the background-subtracted light curves are fitted with equation (22) of Paper 
I, and then we normalize and re-scale the data of the background-subtracted light curves, 
using the method adopted above, with the corresponding fitting curves (the magnitude and 
FWHM position of the background-subtracted light curves data could be well estimated 
from these fitting curves). We find that all pulses in our sample exhibit a reverse S-feature 
deviation from the marginal decay curve, where the (central) profiles of the light curves are 
indeed well below the marginal decay curve beyond the FWHM position of the decay phase 
(see Figs. 4 and 5). When compared with the standard decay form, all of them except two, 
# 3257 and # 5495, show the reverse S-feature deviation as well (for the details see Fig. 5). 
A list of the deviation areas, A\ and A%, from the marginal decay curve and the standard 
decay form, are shown in Table 1, and their distributions are presented in Fig. 6. 

Comparing Fig. 6 with Fig. 3 we find that for most pulses of the sample the widths of 
the decay phase of their corresponding local pulses are sufficiently large (larger than 0.1) 
that they are no more sensitive to the two areas A\ and A^. This suggests that, compared 
to the curvature delay time of the emitting region of the shock (see what discussed in last 
section), the intrinsic pulse decay times are long. This is because the FRED shape due to a 



local ^-function pulse has a characteristic duration set by the curvature time delay — the 
time delay between the arrival times of two simultaneously emitted photons, one from the 
line-of-sight and one from 9 max , but that due to a local pulse with a sufficiently long decay 
time has a different characteristic. At a certain observation time, only photons emitted from 
a limited area could reach the observer in the case of the local 5-function pulse (in this case, 
the area could be marked by 9 — 9 + A0, when A9 is extremely small), while in the case 
of the long decay time local pulse, photons emitted from all the areas concerned could be 
observed (they are emitted from different local times) (note that, in the case of the suddenly 
dimming local pulse, the corresponding area would decline with time). It seems that it is 
this difference that leads to the variance of the light curve characteristic seen in the two 
cases. Under this interpretation, we suspect that the fact that many BATSE pulses are in 
the saturation regime suggests that constraints may be placed on the angular size and radius 
of the emission region (which might deserve a further investigation). 



4 DISCUSSION AND CONCLUSIONS 



A reverse S-feature deviation of profiles of light curves arising from local pulses containing 
a decaying tail from those associated with very narrow or suddenly dimming local pulses, 
the so-called standard decay form, could be seen in Fig. 3 of Paper II. We investigate in this 
paper if FRED pulses observed bear indeed this feature, and if they do how to measure the 
deviation. We define two areas A\ and A 2 to describe the deviations within and beyond the 
FWHM position in the decay phase, respectively. Suggested in our analysis, different from 
the standard decay form, there is a marginal decay curve which reflects the profile of the 
light curve arising from a local S function pulse with a mono-color radiation. We employ a 
sufficiently large sample of FRED pulses of GRBs to study this issue. The study shows that 
the reverse S-feature indeed exists in the profiles of all the pulses concerned when compared 
with the marginal decay curve, while 73 out of 75 individual pulses show the feature as well 
when compared with the standard decay form. We also find that the values of Ai and A 2 
for most pulses of the sample are quite large which suggests that the corresponding local 
pulses might contain a long decay time relative to the time scale of the curvature effect. 

For the two exclusive events, # 3257 and # 5495, the deviation from the standard decay 
form is "positive" , rather than "negative" , beyond the FWHM position in the decay phase 
(say, t' > 0.2). We do not know what causes this exception. Here are several outlets we can 
think of. The first is associated with the background subtracting. We find that over or less 
subtracting the background count would shift the corresponding profile under or over the 
standard decay form beyond the FWHM position in the decay phase. Illustrated in the left 
panel of Fig. 7 is this effect which is quite significant. The second is the impact of the rest 
frame radiation form. Note that the standard decay form is defined when the indexes of the 
rest frame spectrum are taken as a = —1 and (3 = —2.25. As shown in Preece et al. (2000), 
the indexes could take many other values. We wonder if different values of the indexes could 
lead to a much different deviation. Shown in the right panel of Fig. 7 is this effect, where two 
modified standard decay forms, for which the rest frame spectral indexes of the standard 
decay form are replaced with others, are presented. It indicates that different values of the 
indexes could indeed change the profile beyond the FWHM position in the decay phase, 
but the deviation is relatively small (compare the two panels of the figure). 

As suggested in Paper II, equation (1) holds in the case of uniform jets. When # max 
is very small (say, # max ~ 1/T), there will be a turnover feature in the decay tail of the 
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Figure 5. -Plots of the fitting curves (the dot line), the standard decay form (the dash line) and the marginal decay curve 
(the solid line) for the 75 individual pulses of the sample. 
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Figure 6. Distributions of A\ and A2 of all the pulses in our sample. The two lower panels are associated with the marginal 
decay curve and the two upper panels correspond to the standard decay form (the two exclusive events are not included). 



light curve (see Paper II Fig. 2). One could check that, when # max is sufficiently large, 
the turnover feature would not be detectable. In this case, a FRED pulse would also be 
observed. Therefore, the conclusion above holds as well in the case of uniform jets, as long 
as the opening angle of jets is not extremely small. 

Hinted by our analysis and the previous works (see Ryde and Petrosian 2002; Paper II), 
one can conclude that the curvature effect would lead to FRED pulses, and for the pulses 
caused by the curvature effect their profiles would exhibit a reverse S-feature deviation from 
the marginal decay curve. Thus we propose to take the marginal decay curve as a criteria to 
check if an observed pulse could be taken as a candidate suffered from the curvature effect. 

An interesting question arises accordingly, which is that, for those non-FRED pulses of 
GRBs, what one could expect. We suspect that, the reverse S-feature might no longer be ob- 
served in these cases and such pulses might be associated with structure jets or the scattering 
ejecta. This might deserve a further investigation (see Lu and Qin 2005 in preparation). 
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made many helpful suggestions. This work was supported by the Special Funds for Major 
State Basic Research Projects ("973") and National Natural Science Foundation of China 
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Figure 7. The left panel is the plot of the effect of the background subtracting, on the deviation from the standard decay form, 
where the solid line is the standard decay form C(t), the dash line represents the curve of C"(r) = C(t) + 0.05C P (which is 
associated with the case of less subtracting the background count), and the dot line stands for the curve of C"(r) = C(r) — 0.05C P 
(which corresponds to the case of over subtracting the background count). The right panel is the plot showing the effect of the 
rest frame radiation form on the deviation from the standard decay form, where the solid line is the standard decay form (for 
which, «o = —1 and (3q = —2.25 are adopted), the dash line is the modified standard decay form when replacing /3q = —2.25 
with Po = —1, and the dot line represents another modified standard decay form by replacing ojo = — 1 with c?o = —0.5 and 
replacing /3q = —2.25 with /3q = —4.5. 



Table 1. A list of the estimated values of the deviation areas, A\ and A2, from the marginal decay curve and the standard 
decay form, for the pulses of our sample 
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Note: Ai (M) and A2 (M) are the two deviations of the profile of a light curve from the marginal decay curve. 
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